BNL-HET-08/9, CERN-PH-TH-2008-081, VPI-IPNAS-08-09 

Feasibility Study for Measuring Geomagnetic Conversion of Solar 
Axions to X-rays in Low Earth Orbits 

Hooman DavoudiaslQ 

Department of Physics, Brookhaven National Laboratory, Upton, NY 11973, USA 

Patrick Huber[| 
Physics Department, Theory Division, 
CERN, 1211 Geneva 23, Switzerland and 
Department of Physics, Virginia Tech, Blackshurg, VA 24062, USA 
(Dated: March 20, 2009) 

We present a detailed computation of the expected rate for Geomagnetic Conver- 
sion of Solar Axions to X-rays (GECOSAX) along the orbit of an x-ray satellite. We 
use realistic satellite orbits and propagation in time. A realistic model for the Earth's 
magnetic field, which properly accounts for its spatial non-uniformity, is used. We 
also account for the effect of the Earth's atmosphere on the propagation of x-rays in 
our calculation of axion-photon conversion probability. To estimate possible sensi- 
tivities to the axion-photon coupling gay, we use an actual measurement of the ex- 
pected backgrounds by the SUZAKU satellite. Assuming a detector area of 10^ cm^ 
and about 10^ s of data, we show that a 2 cr limit of ga-y < (4.7 — 6.6) x 10^^^ GeV^^ 
from GECOSAX is achievable, for axion masses rua < 10~^ eV. This significantly 
exceeds current laboratory sensitivities to ga-y 



I. INTRODUCTION 

Weakly interacting light pseudo-scalars are well-motivated in particle physics. For ex- 
ample, experimental observations require the size of CP violation in strong interactions, 
parametrized by the angle 6, to be quite small: 6 ^ 10^^°. However, the symmetries of the 
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Standard Model (SM) allow ^ ~ 1; this is the strong CP problem. An elegant solution to 
this puzzle was proposed by Peccei and Quinn (PQ) 0, where a new U{1) symmetry, 
anomalous under strong interactions, was proposed. This U{1) symmetry is assumed to be 
spontaneously broken at a scale fa, resulting in a pseudo-scalar Goldstone boson a 
the axion. Non-perturbative QCD interactions at the scale Aqcd ~ 100 MeV generate a 
potential and hence a mass mf*^ ~ Aqq^Z/q for the axion. Experimental and observational 
bounds have pushed fa to scales of order 10^ GeV or more. As fa sets the inverse coupling of 
the axion to the SM fields, the current data suggests that axions are basically 'invisible' and 
very light. Here we note that some cosmological considerations related to the over closure of 
the universe suggest an upper bound fa ^ 10^^ GeV for the PQ axion |^, 0, §. Apart from 
the considerations related to the strong CP problem, axion-like particles are ubiquitous in 
string theory. In addition, axion-like particles have been used in various astrophysical and 
cosmological models^. In the following, the term axion is generically used to refer to any of 
the above, or other, possible instances of such weakly interacting hght pseudo-scalars. 
The coupling of the axion to photons is given by 

jCay = = ga^aE-B, (1) 

where M ~ {7i/a)fa and a ~ 1/137 is the fine structure constant. F^^, is the electromagnetic 
field strength tensor, F'^'^ is its dual, ga^y = M-^ is the axion-photon coupling; E and B 
are the electric and magnetic fields, respectively, corresponding to F^^^. The interaction in 
(H) makes it possible for hot plasmas, like the Sun, to emit a flux of axions through the 
Primakoff process ||lT|. This same interaction has also led to experimental proposals ||12| for 



detecting the axion through its conversion to photons in external magnetic fields. Various 
experimental bounds, most recent of which is set by the CAST experiment [|13|, suggest that 
da-r ~ 10^^° GeV^^. For a review of different bounds on axion couplings, see Ref. [Q. 



In what follows, we study the feasibility of a recently proposed approach for detecting 
solar axions with an x-ray telescope in orbit |T^, based on geomagnetic conversion of solar 
axions to x-rays (GECOSAX)^. The estimate of the expected x-ray flux presented in Ref. |15 
was based on a number of simplifying assumptions: 



^ See, for example, Ref. 1^. 

^ The possibility of using planetary magnetic fields as a conversion region for high energy cosmic axions 
was discussed in Ref. [[16|. 
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1. 



The satellite orbit was a circle. 



2. The orbit was aligned to lie in the equatorial plane of the Earth. 

3. The Earth axis was perpendicular to the Ecliptic. 

4. The available conversion length was taken to be the altitude of the satellite. 

5. The magnetic field was assumed uniform and perpendicular to the direction of axion 
propagation. 

6. The effect of the Earth atmosphere was neglected. Consequently, the effective mass of 
the photon in medium was ignored; — » 0. 

These assumptions allowed a treatment of GECOSAX within the same formalism relevant 
for helioscopes |12, |T^, i.e. axion-photon conversion in vacuo in a constant magnetic field 



which is perpendicular to the direction of the axion momentum. The conversion rate , in 
the limit of vanishing axion mass rria ^ 0, is then simply given by 



Taking gay = 10-1° GeV-i and using CAST parameters B = 9T and L = 10 m, we obtain 
P^^ ^ 2 X 10^1''. Replacing the CAST magnet with the geomagnetic field and taking L to 
characterize a low-Earth-orbit, we have B = 3 x lO^^T and L = 6 x 10^ m and thus get 
for P^^ ^ 8 X 10^1^, which is only a factor of about 25 smaller than the CAST conversion 
probability. However, in Ref. |]l5l, it was noted that this can typically be overcompensated 
by the larger detection area of an orbiting x-ray telescope. In |jl5| it was shown that the 



resulting x-ray signal on the dark side of the Earth would have a number of unique features 
which would make it very hard to be mistaken for anything else: upward going x-rays, 
T = 4keV black body spectrum, direction from within 3' from the center of the Sun and 
characteristic modulation with B^LF' . 

In the following, we would like to address the simplifying assumptions one by one and 
illustrate their effect on the actual axion-photon conversion rate. In section y, we will 
discuss how to account for the proper satellite- Sun-Earth geometry. We will also discuss 
how to compute the actual position of any given satellite. This will address assumptions 
1-4. Next we will address assumption 5 in section |T| with a proper magnetic model. In 
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section we will present a full treatment of axion-photon conversion in a dispersive and 
absorptive medium, using a model of the Earth's atmosphere, which addresses assumption 6. 
In section |V], we will present the resulting x-ray fluxes for various satellites, followed by a 
discussion of achievable sensitivities for Qa-y, in section |V|. Finally, we will present a discussion 
of our results and the future outlook, in section |V14 

II. GEOMETRY AND SATELLITE ORBITS 

The basis for any detailed calculation of the axion-photon conversion rate is a correct 
description of the geometry. We introduce the following notation: plain capital Latin or 
Greek letters denote a point, where O denotes the origin of our coordinate system. If P is 
a point, then its position vector OP is denoted by p. For any vector v its length is denoted 
by V. The unit vector along the direction of v is denoted hy = v/v. The Sun is at S, 
the center of the Earth is at E and the x-ray satellite is at X. For any point P, we define 
its height vector hp as the vector which goes from the surface of the Earth to P and is 
perpendicular to the tangent plane of the Earth surface at its starting point. The starting 
point is called the footprint of P and denoted by Pp. This definition may seem involved, 
however, it also holds for the actual geoid and not only for a spherical Earth. 

The abstract definition of the problem can be easily done without specifying a coordinate 
system. The actual numerical calculation, of course, has to specify a definite coordinate 
system, which is relegated to appendix 0. 

In Fig. |l] a two dimensional schematic view of the problem is shown at one instant in time 
to- Since axions/photons travel close to/at the speed of light c, they experience an essentially 
static environment. The geometry only changes notably on timescales large compared to the 
propagation time Tp ~ 10^ m/c ~ 10"'^ s. Therefore, the geometry can be regarded as fixed at 
any given time to- The first task is to determine the path traveled by the axions. The axions 
propagate in a straight line from the Sun 5* to the satellite X. For axion conversion, however, 
only the part of SX which is on the dark side of the Earth is relevant. Thus, the intersections 
of SX with the surface of the Earth have to be found. We account for the ellipsoidal shape 
of the Earth and use the so called WGS72 parameters^. In general there can be none, one or 

^ For polar orbits the difference in radius of a spherical and ellipsoidal Earth can cause up to a 10 s difference 
in the duration of the dark orbit. 
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FIG. 1: Geometry of the GECOSAX configuration drawn in the plane spanned by the center of 
the Earth, the center of the Sun and the sateUite's position at to. 

two such intersections. The solution on the dark side will be denoted by L. This allows us 
to define the line of sight (LOS) a = LX. Note, that this definition of being on the dark side 
is purely geometrical and neglects the angular diameter of the Sun, atmospheric refraction 
and absorption. The angular diameter of the Sun reduces the useful part of the orbit by 
about 4 s. This follows from the fact that the Sun subtends 0.5° and the sateUite travels 
a full circle in about 90 minutes. The effect of the Earth atmosphere is quite a bit larger 
since it becomes non-transparent for x-rays below an altitude of about 50 — 100 km, thus 
increasing the effective radius by that amount, which by explicit calculation would increase 
the dark orbit duration by a few times 10 s. This ovcrcompensates for our neglect of the 
solar diameter and makes our overall treatment conservative. Using the position of L as 
defined above, one can parametrize the position of any point along the line of sight 

p(A) = Ae; + r. (3) 

In the course of the actual calculation the height hp oi P is needed, since the air density is 
a function of the actual height. Note that the height of the satellite hx is always smaller than 
the length of the LOS, i.e. hx < ce. This implies that assumption 4 is in fact conservative, 
and we will find that the relative increase of a with respect to hx will compensate largely for 
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the losses in x-ray flux due to the other effects considered in the following. The algorithm 



for the solar position is taken from |T8[ and it's accuracy is better than 1'. 

The idealized orbit of any satellite is a solution to the Kepler problem. Thus, knowing 
the satellite's position and velocity at time t^, it is possible to predict its future position at 
ti- In reality, there are, however, various factors which lead to deviations from the simple 
Kepler orbit, among which are: non- vanishing higher multipoles of the mass distribution 
of the Earth, atmospheric drag, gravitational influence from the Moon (and to a lesser 
degree from the Sun), etc.. The prediction of satellite orbits is, of course, a matter of 
great importance for operators of satellites and is also needed for military applications. The 
aforementioned effects disturbing the simple Kepler orbit can be accounted for in a general^ 
perturbation theory. Many of the perturbations are well know and can be quite exactly 
computed. From the observation of the actual position and velocity of a satellite at time to 
it is possible to extract the unperturbed Kepler orbit, which would follow in the absence of 
any perturbing factors. In order to obtain an accurate prediction for the future, a specific set 
of perturbations is taken into account. In doing this, it is crucial that the initial unperturbed 
Keplerian orbit data is extracted using a model which is compatible with the algorithm used 
for future positions. 

One standard format is the so called 'NORAD^ element sets', and a description of the 
perturbation model called SGP4^ can be found in [IT3. Since satellite propagation is done 



in perturbation theory, errors inevitably will accumulate and render the predictions unre- 
liable. Therefore, element sets for basically all active satellites are issued periodically by 
NORAD and made accessible at ||2^ in the so called 'two line element' (TLE) format. The 
implementation of the NORAD orbit prediction algorithm we use is taken from the predict 



program, which is an open-source C language satellite tracking software [21] . It directly 



General in the sense, that the resulting theory is applicable to a wide, general class of orbits and not 
restricted to particular orbits like e.g. ones with a low eccentricity. 
^ NORAD is the North American Aerospace Defense Command. 

^ SGP4 stands for 'simplified general perturbation version 4'. Historically, one distinguished SGP4 and 
SDP4, where the latter one is used for 'deep-space' orbits with periods longer than 225 minutes. For most 
parts of this work, we use only SGP4, i.e. orbits with periods smaller than 225 minutes. 

^ In reality, it seems that all implementations found in open accessible sources go back to various, different 
original implementations by T.S Kelso pO| | in a number of programming languages. The one we are using 
is no exception. 
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takes the TLE of a satellite, a time to and returns its position in the ECI^ at to. All satellites 
are indexed by NORAD using so called US SPACECOM identification numbers, these are 
5 digit numbers starting with 00001 for the SPUTNIK satellite. We will use these 5 digit 
numbers to refer to all satellites in this paper, any satellite names are written in capitals. A 
list if all US SPACECOM IDs and the corresponding names is given in table |ITT[ 

Some remarks about our use of TLEs and SGP4 are in order. We use SGP4 since it 
is the simplest general purpose algorithm and the necessary input data, the TLEs, are 
easily available. SGP4 is by no means state of the art, it was developed to allow reliable 
tracking of thousands of objects with the limited computing power available in the 60s and 
70s. Clearly, in an actual experiment one would use telemetry data and direct numerical 



integration, possibly even GPS, thus reducing any position errors to around 100 m or less p2 
The accuracy of predictions made with SGP4 relative to GPS position determination was 
studied in detail in p2|. The typical errors for SGP4 are about ±5 km cross-track, i.e 



perpendicular to the satellite momentum, within ±15 days from the epoch of the used TLE, 
whereas the in-track error, i.e. along the orbit, is about ±20 km. For a satellite moving 
at about lOkm/s this gives rise to a timing error of about 2 s. We checked that these 
errors have, in fact, a very small impact on the average axion signal, since the satellite still 
reaches every point under more or less the same circumstance with respect to magnetic field 
orientation and direction to the Sun. For this test, we used historic TLEs of satellite 27370 
(RHESSI), obtained from [^, issued about 8 months apart. We found that there was a 
time difference of several minutes in when the satellite entered the dark orbit, but once we 



corrected for this time shift, the GECOSAX fluxes where identical to within ~ 10%. In ||22 
it is shown that the inter-TLE variation is a good indicator for the actual accuracy. This 
is clearly an extreme example since TLEs are re-issued about every other week. Thus, we 
conclude that SGP4 with current TLEs is accurate enough by a large margin. 

To summarize, we found that the corrections due to the proper treatment of the geometry 
are quite large and lead to a pronounced variation of the expected flux along each orbit 
because of changes in the length a of the LOS. Geometry related effects are accounted 
within about ±25 km or ±5 s in our calculation, which introduces less than 10% error in 
the GECOSAX signal prediction. The errors introduced by our simplified treatment of 



This is an Earth Centered Inertial set of coordinates, discussed in appendix ^ 
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FIG. 2: Map of the total magnetic field strength at sea level for 2008.5 [^]. Red, thick contours 
are in steps of lOOOOnT, thin red ones are in steps of 2 500nT. The number on top of the red, 
thick contours are the magnetic field strength in nT. The black dots denote the positions of the 
magnetic dip pole for each year from 2005 till 2010. The coordinates give the position of the dip 
pole in 2008.5. The map is a Winkel tripel projection. 

coordinate transformation in appendix ^ are negligible in comparison to the intrinsic errors 
of the satellite orbit prediction. 



III. MAGNETIC FIELD OF THE EARTH 



To a first approximation, the magnetic field of the Earth is a dipole whose axis intersects 
the surface at the magnetic poles, which do not coincide with the geographic ones. This 
mis-alignment of rotation and magnetic axes alone would induce a typical periodic variation 
of the x-ray flux produced by GECOSAX. However, the geomagnetic field has various other 
irregularities and deviations from a simple dipole form. 

We use a realistic 3-d model of the magnetic field of the earth, the so called World 
Magnetic Model 2005 p4|, which is available in machine readable form at p3|. The World 



8 



Magnetic Model is the standard model of the US Department of Defense, the UK Ministry 
of Defense, the North Atlantic Treaty Organization (NATO), and the World Hydrographic 
Office (WHO) for navigation and attitude/heading referencing systems. It is intended to 
meet its stringent error specifications (better than 1%) from sea level up to an altitude 
of 600 km. Since it is given as a series expansion in spherical harmonics, mathematically 
it stays well-defined out to larger radii. If there were no electrical currents in the upper 
atmosphere and no solar wind, i.e. additional sources of magnetic fields outside the Earth, 
the model would be accurate up to many Earth radii. In practice, interactions with the solar 
wind and atmospheric electrical currents produce magnetic fields of around 100 — 500 nT, 
at an altitude of ~ 1000 km, during magnetically quiet times and perturbations can reach 
up 2000 nT during strong magnetic storms. Magnetic perturbations are indexed by the 
Ap index, which is the daily average of the ap index. It denotes the deviation from the 
most disturbed component of the local magnetic field vector from its mean, undisturbed 
value in units of 2nT. The Ap index is derived from the observations of 11 geomagnetic 
observatories and has been regularly collected since 1932. Ap values larger than 100 are 
classified as indicating a severe magnetic storm. Only 1% of the days from 1932 till 1992 
have reported a value of Ap > 100 Hence, for almost all observation conditions the 

errors introduced by the day to day variability of the geomagnetic field will be small. Thus, 
the errors introduced by using the magnetic model up to an altitude of 1000 km are certainly 
less than 10%, most likely much less than 5% Therefore, in principle, it seems 

feasible to extend the permissible range of altitudes maybe up to 1-2 Earth radii, however, 
this would require a more careful analysis of the external magnetic fields, which is beyond 
the scope of this work. We, therefore, will restrict all analysis to altitudes below 1 000 km, 
unless otherwise mentioned. 

The magnetic model also includes a prediction of the annual variation of geomagnetic 
parameters from 2005 to 2010. From these variations we expect a less than 1% annual 
change in the relevant parameters. Given this 3-d vectorial map, we compute the transverse 
S-field along the axion path. The total field strength is shown in Fig. 
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IV. AXION PROPAGATION 



The axions and x-rays will have to traverse the upper Earth atmosphere, which causes 
absorption and refraction of x-rays and hence will also influence the axion-photon conversion 
probability. To a rough approximation the interaction of x-rays with an energy of few 
keV with air can be described by Thomson scattering from free electrons |2^. Air mostly 



consists of nitrogen and oxygen having atomic numbers Z of 7 and 8, respectively. The 
binding energies of the innermost electrons thus are about 13.6 eV ~ 600 — 800 eV and 
thus small compared to the photon energy for most of the range of interest. On the other 
hand, the photon energies are very small compared to the rest mass of the electron and 
hence the scattering is highly non-relativistic and pair creation cannot take place. Thus, all 
effects on x-ray propagation should be a function of the electron density, which itself closely 
traces the mass density of air. At standard temperature and pressure (STP) of 273.15 K 
and 101 325 Pa, we use as volume (molar) fractions 78.1% N2, 21.0% O2, and 0.93% Ar^. 

The absorption length A = F"^ for x-rays of energy 1 — 10 keV, in air at STP, has been 
obtained from ||2^, 0. We have assumed that x-ray absorption scales with the electron 
density along the axion path. Assuming a constant composition of the atmosphere with 
altitude, the electron density^^ is directly proportional to the mass density. 

To include refraction, we use the effective photon mass given by |1^ 



= 47rro[p/i/(Am„)], (4) 

where vq ~ 2.82 x 10^^^ m is the classical electron radius, A is the atomic mass number of 
the gas (atmosphere), rriu is the atomic mass unit, p is the gas density, /i ~ Z, and Z is 
atomic number of the gas. This formula can be generalized for a compound gas, like the air, 
by noting that the quantity pZ/ {Arriu) is the electron number density Ue for the medium, 
which in the above equation is assumed to be made up of only one element. 

For a simplified derivation we assume that air is composed of 78% N2 and 22% O2, 
by volume. In the ideal gas limit, the volume and molar fractions are the same. Given 
that for dry air the molar density is pair — 44.48 mol m~'^ (STP), we find that Ug — 



^ With Z = 18, Ar has binding energies in the keV range. However, due to its small molar fraction we may 
ignore it for the exposition, here. In the numerical analysis it is accounted for. 
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We will comment on this in more detail, later in this section. 
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44.48 A^a(0.78 x 14 + 0.22 x 16) =i 3.87 x 10^*' m-^ where Na is Avogadro's number. Then, 
Eq. (ID yields 



m. 



\ 1/2 



^-°-^Mk^) (5) 



We checked that using the full energy dependence of /i given in |g9|, |32] and a more detailed 
air composition does not change by more than 2%. 

As the following general argument will show, the axion conversion path length is only 
logarithmically sensitive to density variations. In all density models, the density profile can 
be locally described by an exponential with an altitude dependent scale height H{h) 

p{h) = p^e-Tm , (6) 



where both po and HiK) can be time dependent. In order to estimate the impact on axion 
propagation, we need to understand up to what altitude absorption plays a role and what 
impact a finite photon mass has. In order to asses the sensitivity to changes in absorption 
we can compute the escape probability pesc of a x-ray photon from a given altitude ri to 
infinity: 

/■oo 

Pesc = e"''*' with Cp = / dxV{x), (7) 

where r(a;) is the inverse absorption length and Cp the so called absorption coefficient. r(a;) 
itself is a function of the density and is given by 

r(x) = pp{x) , (8) 

where p is the mass attenuation coefficient. We now can define the escape altitude a^sc by 
demanding pesc = 1/e. This then translates to the following condition 

l = Cp = / dxT{x) = / dxppoe"^. (9) 



CSC 



Assuming that H{x) = H = const, we can easily solve for x and obtain 

7]esc = Hln{Hpop). (10) 

The effect of a non-vanishing is the same as the one of a non-vanishing axion mass: 
they both affect the oscillation or coherence length for photon-axion conversion. The con- 
version probability will be suppressed whenever the oscillation length is short compared to 
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the available path length. For the following discussion we define the oscillation length Ij to 
be due to the finite photon mass 

= (11) 

where lo denotes the photon energy. In order for to be negligible, we require < L 
1000 km. Since oc y^p, it follows that /-y oc p~^. The altitude rj at which reaches a 
certain value can be computed, using Eqs. |^ and || 



7] = Hhi 



1.65 X lOM 7-^ 



(12) 



^ km J \kgm ^ J \ iu 

Let 7]^ denote the value of rj at 1^ = 1000 km. Again we find that the dependence of rj^ on 
Po is only logarithmic. Therefore, all atmospheric effects will depend only weakly on the 
precise value of density. 

The Earth atmosphere is roughly compromised of three layers: the homosphere from 
— 90 km, the thermosphere starting at 90 km and ranging up to 250 — 400 km depending 
on solar and geomagnetic activity, as well as the exosphere which begins at the top of 
thermosphere and extends into space. In the homosphere, winds and turbulence mix all 
species very well and thus the composition is independent of altitude. In the thermosphere 
turbulent mixing ceases to be effective and the different species start to diffuse separately. 
This diffusion is driven by gravitation and thermal gradients. In the exosphere, finally, 
the mean free path of the lighter atoms like hydrogen becomes large enough such that 
they can escape to space. In the thermosphere solar energy is absorbed by the photo- 
dissociation of molecular oxygen. Thus, a sizable amount of free, atomic oxygen appears 
and our assumption of constant composition fails. However, in terms of electron density, 2 
oxygen atoms have the same number of electrons as 1 oxygen molecule, therefore the total 
mass density is still a very precise indicator of the electron density^^. 

There are no weather phenomena, in the ordinary sense, in the thermosphere, nonetheless 
its density does depend on various variable factors. This density distribution depends on the 
energy input from the Sun via extreme ultraviolet light (EUV) and due to direct heating by 
charged particles from cosmic radiation and solar wind. As a result, the density depends on 
the amount of solar EUV radiation, which itself depends on geographic latitude, the time of 
the year and the apparent local solar time. Also the solar activity and geomagnetic activity 
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We verified that this holds to better than 1% in the relevant altitude range. 
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FIG. 3: Total mass density p as function of altitude. The black line shows the default density 
profile used throughout this paper and corresponds to the medium solar activity case recommended 



in [28|. The green/gray band denotes the maximal excursions from this default density profile 
predicted by the NRLMSISE-00 atmospheric model . The inset shows the photon mass induced 
minimal oscillation length (blue/dark band) and the absorption coefficient Cp (green/light band) 
as functions of the altitude in the relevant range. The bands are due to full variation of density 
profiles shown in the main plot. The arrows at the bottom denote the main cause for the density 
variation in that altitude range. 

as well as the location within the geomagnetic field do have a non-negligible influence. For a 



general overview see |28|, Direct measurements of the total mass density in the altitude 
range from 150 — 200 km can be performed by the observation of the decay rate of very low 
altitude satellite orbits. These measurements indicate densities of a few times 10~^ kg m~^ 
at 150 km. The observed diurnal variations are 25%, whereas the observed seasonal variation 



is somewhat larger with 40% ^| 
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In the MSISE-90 ||3^ model^^ is recommended for use in space missions, with three 



average density profiles from — 900 km corresponding to three different levels of solar and 
geomagnetic activity. We use the one called medium activity, corresponding to -F10.7 = 
F^Qj = 140 and Ap = 15. Where, -F10.7 denotes the A = 10.7 cm flux density from the Sun 
in units of 10^^^ Jm~^. -Fff^f is the 81 day average of -Fio.7- The -F10.7 index closely traces 
the solar UV emissions and the sun spot number. Ap is the daily average of the ap index, 



which measures the perturbation of the geomagnetic field in units of 2nT, see section |T1. 
The corresponding density profile is shown as black line in figure ^. 

In order to study the impact of variations of the density profile we use an updated 



version of MSISE-90, the so called NRLMSISE-00 |31| model. The differences between the 



two models for the same set of input parameters, are however small. NRLMSISE-00 takes as 
input the day of year, the local apparent solar time, the geodetic latitude and longitude, -F10.7 
and -Fft^f and Ap (or a series of average values of Ap). In order to to estimate the maximal 
possible density excursions, we varied: -F10.7 and FfJ^f jointly from 40 to 380^^, Ap from 
to 100^^, longitude from 0° to 90°^^, the day of the year from 1 to 365, the local apparent 
solar time from h — 24 h. For each altitude we determined the minimal and maximal value 
of density due to all these different input parameters, the result is shown as the green/gray 
band in figure ^. In the altitude range from 50 — 135 km seasonal changes and the geodetic 
latitude have the greatest effect, wheres for higher altitude the main effects are due solar and 
geomagnetic activity. These regions are indicated by black arrows in figure ^. The obtained 



values for density variation agree well with the ones found in [^, The inset in figure 



shows how T^esc and r]^ change due to those density variations. The values and ranges are: 

^csc = 78^*^3 km and rj^ = 118 ± 3 km. (13) 

We see that the limiting factor is indeed refraction and not absorption and both factors need 
to be included for an accurate calculation. Assuming an axion conversion path of around 
1000 km, this is less than a 1% change. We actually verified that the GECOSAX flux does 



MSISE stands for 'mass spectrometer and incoherent scatter, extended'. 
^•^ This covers the extremes during one full 11-year solar cycle, according to [ p8| . 

In principle, Ap can reach values of up to 300 during the strongest geomagnetic storms. These times, 
would however have to be discarded anyway since the fidelity of geomagnetic model can not be ensured at 
these times. Note, that the density variations caused by Ap > 100 especially affect the vicinity of 120 km, 
i.e. precisely the regions where 77^ will be located. 
The —90° to 0° range just swaps result between summer and winter. 
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not change by more than 5% due atmospheric density variations, therefore the inclusion of 
atmospheric effects via an average density profile is fully warranted. 



A. Axion conversion probability 

The probability for axion-photon conversion including the full path and medium depen- 
dence is given by fl^ 

2 
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and 



A = ^ exp |- 1^"' dX r [/ip,(A)]} (15) 



BU\',t)^\B[pt{\')]x4\. (16) 

Here, uj is the axion energy. The time dependence of Pa-y is entirely due to change of the 
geometry with time as explained in section For each time t we solve for the position of 
the satellite Xt and for the position of the Sun St- This information is used to derive Pt(A), 
the parametric form of the axion path or the line of sight as defined in Eq. |[ The quantity A 
parametrizes the position along the line of sight and at denotes the length of the line of sight 
for the time t; we have A G [0, at]. F and only depend on the density of the atmosphere 
which itself is a functions of the height above mean sea level hp^{X). Only B has a complete 
dependence on the position vector p. For the various necessary coordinate transformations 
we refer the reader to appendix 0. 

The integral in Eq. ( p^ has no closed form solution and therefore has to be integrated 
numerically. For the numerical integration a problem arises at very low altitudes, where 
the air density is high and hence is large. This gives rise to extremely fast oscillation 
of the integrand of the innermost integral, i.e. the integration with respect to A" becomes 
practically impossible for sufficiently small heights. On the other hand, F also becomes 
very large and thus those parts of the path where is very large do not contribute to the 
transition amplitude. The solution is to reverse the direction of integration by substituting 
A with at — X and at the same time all integrals now run from at to 0. Next, the integral 
is partitioned using a simple bisection: First the integral from at to at/2 and then from 
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ttj /2 to at /4 . . . , until the contribution of the last part evaluated is smaller than a preset 
precision goal, in our case this is AP/P = 10~^. 



B. Solar axion flux 



The Sun produces axions throughout its whole interior, although the hottest regions with 
the highest photon density contribute the bulk of the axion production. Since the angular 
size of the Sun and the axion producing region are non-negligible compared to the typical 
angular resolution of x-ray telescopes, we cannot treat the Sun as point source of axions. 



A detailed study of solar surface axion luminosity has been presented in |T3|. Its results 



have been made accessible to us in machine readable format by one of the authors [36 



We will denote the solar surface axion luminosity by (pa{r, E), following the notation in [IT3 



where r is the dimensionless fraction of the solar radius. Thus, the solar axion flux spectrum 
produced by the Sun up to a certain radius is obtained by 



■0 



dE 
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drr ipa{i^, E) . 



(17) 



The usually quoted solar axion flux assumes = 1 and can be written as |13[| 



d^a 







dE 



6.02 • 1010^2-481 . 



10-10 GeV 



-1 



cm 



^^keV 



-1 



(18) 



Throughout this work, unless otherwise stated, we use ga-y = 10-10 GeV-i. Note that the 
vast majority of this flux originates within the inner 20% of the solar radius, i.e. Vg = 0.2. 
Since the background will be proportional to r^, the signal significance will not be optimal 
for Ts = 1, but for some smaller value. This issue is studied in detail in appendix and we 
adopt Tg = 0.13. The flux we are using is then given by 











dE 



dE 



1.72- 1010^3-210. e-^A-135 



T-.=0.13 



10-10 GeV" 



cm 2 s ikeV i 



(19) 



The resulting axion fiuence in the energy range from 1 — 10 keV is 3.58 ■ lOH axions cm-2 3-1 
for = 1 and it is 2.21 ■ lOH axions cm-2 3-1 _ q Also, the mean axion energy 

changes from 4.2 keV for = 1 to 4.8 keV for = 0.13. In appendix it is demonstrated 
that the loss of about 40% in signal is compensated by a large decrease in background. 

From a comparison of the results obtained in and in [13 1, it is estimated that this flux 
calculation is accurate to about 5%, due to changes in the underlying solar model. In our 
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analysis, we will also neglect the annual variation of the Sun-Earth distance which causes a 
5% modulation of the signal. 



We now have all the tools at hand to study the GECOSAX effect in detail for any given 
satellite. Before we delve into finding the optimal orbits in the following section, we briefly 
describe the GECOSAX flux along a typical orbit. This will help to clarify some notation 
and to give a basic overview of the issues involved. Figure |^ shows one dark orbit, i.e. that 
part of the orbit which is in the Earth shadow, of satellite 25399 beginning at 2007-12-31, 
23:48:33 UTC and lasting 1524 s. We will call the beginning of a dark orbit and the end 
t^. The duration of the dark orbit is then given by — t^. Since the satellite may not 
be able to start observation immediately at t^, we may have to cut away some time at the 
beginning and end of the dark orbit, called tcut; thus the useful duration of the orbit is 
tu = t{ — t]^ — 2tcut- Generally, the duration of dark orbits will vary throughout the year and 
is different from orbit to orbit. For some parts of the year, there may even be no dark orbits 
at all, quite analogous to the fact that during summer the Sun does not set north of the 
polar circle. The integrated useful GECOSAX flux for each dark orbit o, which corresponds 
to the blue shaded area in panel (a) of Fig. ^, is given by 



where t^J and Pa-^it'^E') are different from orbit to orbit and have to be computed corre- 
spondingly. 

Panel (b) of Fig. ^ shows that the axion conversion path is always considerably longer 
than the altitude of the satellite, its length ranges from about 1.6 times up to 4 times the 
altitude of the satellite. This fact is responsible for most of the signal increase with respect 
to our previous estimate in [jl5| . Also, within the first and last few 10 s of the dark orbit 
there is large variation in this length: nearly one half of the total path length variation 
happens within the first and last 60 s of the dark orbit. Therefore, in order to be not overly 
sensitive to errors in timing and positioning we exclude those first and last 60 s from the 



V. RESULTING X-RAY FLUXES 



A. 



General orbits 




(20) 
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600. 800. 
t-td [s] 



FIG. 4: A typical orbit for satellite 25399. tj; is 2007-12-31, 23:48:33 UTC and the duration of 
the dark orbit is 1524 s. Panel (a) shows the resulting x-ray flux from GECOSAX integrated over 
the energy range from 1 — lOkeV. The blue shaded area gives the useful integrated GECOSAX 
flux after accounting for tent- Panel (b) shows the length of the axion conversion path in units of 
the actual altitude of the satellite (line labeled LOS) and the total magnetic field strength at the 
location of the satellite B in units of 10 ;uT (line labeled B). The gray shaded areas are the first 
and last 60s of the dark orbit. 



analysis. This is indicated by the gray shaded areas in Fig. |^. We also see that the variation 
of the magnetic field is non- negligible. This will depend strongly on the path of the satellite 
with respect to the geomagnetic field. The orbit of satellite 25399 has an inclination of 98° 
and thus the satellite does traverse the region of the geomagnetic poles; in this case it is 
the south geomagnetic pole. Also, Fig. ^ shows that the orientation of the magnetic field 
is nearly parallel to the axion path at the beginning of the dark orbit since the very large 
path length right at beginning does not cause a corresponding increase in GECOSAX flux. 
At the end of the dark orbit the magnetic field has a larger component perpendicular to the 
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axion path and hence the increase in path length is well reflected in a corresponding increase 
in the GECOSAX flux. 



B. Optimal orbits 

In this section we will apply the formalism developed in the previous sections to determine 
what constitutes an optimal orbit for observing GECOSAX. In some approximation, the sig- 
nal is proportional to B^L^, therefore we would like to have orbits which have the maximum 
possible path length in the highest possible magnetic fleld. This points to high altitude 
satellites traversing the region of geomagnetic poles. This requires inclinations greater than 
about 70°. Instead of designing an optimal orbit, which given the many free parameters is 
a daunting task, we took a sample of existing orbits, i.e. orbits which actually are used or 
have been used for real scientiflc satellite missions. We took the TLEs of 50 satellites with 



apogees below 1000 km from |2^. Nearly, all of these orbits have low eccentricity < 0.1. 
The apogees of these satellites have an approximately GauBian distribution with a mean of 
650 km and a standard deviation of about 120 km. The inclinations are strongly clustered 
around 80°; 28 satellites have inclinations in the range 70° — 90°. Thus, this sample seems 
to be well suited for our study. The US SPACECOM identiflcation number and the number 



of the TLE set used are given in Table |T|. Our goal here is solely to determine the most 
suitable orbit and not the most suitable mission, which is the combination of satellite and 
orbit. Thus, in the following if we speak, say, of satellite 25544^^ we actually just refer to 
its orbit and not the instruments or the satellite itself. 

In determining the optimal orbit we need to distinguish two different observational strate- 
gies, called 'turning mode' and 'flxed mode'. A turning mode satellite needs to avoid direct 
exposure of its x-ray detection system to the sunlight (visible and x-ray) in order to prevent 
any permanent damage. Typically, this Sun avoidance angle is about 30° and maximum sus- 
tained slew rates of 6° min~^ have been demonstrated. Thus such a satellite enters the Earth 
shadow pointing 30° away from the Sun, then it needs 30/6 = 5 min to turn into observation 
position pointing to the Sun. Since these numbers may vary from mission to mission and 
some safety margin will be necessary, we will discard the flrst 10 minutes at entry into the 
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This is the International Space Station (ISS). 
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Earth shadow as well as the last 10 minutes prior to exit of the Earth shadow, thus giving up 
20 minutes of each orbit. Fixed mode satellites have instrumentation which can withstand 
direct irradiation by the Sun, or protective shields that can be deployed quickly (compared 
to slew time), and thus can do their maneuvering to point to the Sun in the bright parts 
of their orbit. Therefore, in principle, fixed mode satellites can use the entire part of their 
orbit in the Earth shadow. We have noted previously that the precise time of entry or exit 
of the Earth shadow are somewhat uncertain due to geometrical and refractive effects. Also 
the axion conversion rate has its peak values at the entry and exit points where its time 
derivative is largest. Combining these factors, the GECOSAX rate within first and last few 
10 s of each dark orbit have fairly high uncertainties. Therefore, we will exclude the first 
and last 60 seconds of each dark orbit from the analysis even for fixed mode satellites. This 
is by no means a technical necessity but is a conservative choice to ensure reliability of our 
results. 

For each of the satellites in Table |T| we computed the integrated, energy averaged 
GECOSAX flux ($)gcx for each orbit (about 5500 per satellite) from January 1st 2008 
till December 31st 2008 with time steps of 60 s^^. ($)gcx is computed in the limit of 0. 
($)gcx accounts for the time cut away at both ends of the dark orbit; as a result ($)gcx is 
different for turning and fixed observation modes. The signal per unit area for a given orbit 
o is then given by cTq = ($)gcx ^ tu- reality there will be some background b as well to the 
measurement and thus the pertinent quantity to optimize is not the total signal s but signal 
divided by the square root of the background s/Vb. We call this quantity significance. The 
signal is given by So = ctq x A, where A is the area of the detector. The background is 
given hj bo = F X A X f^, where F is the background rate per unit area. F here is to be 
understood as the background rate / integrated over x-ray energies from -Emin = 1 keV to 
Emax = 10 keV and over the solid angle covered by the axion producing region in the Sun. 
Assuming a uniform distribution of the background in energy and solid angle, we obtain 



the following relation F = f2s(ii^max — E^i^f. As mentioned in section [IV B| and explained 
in detail in appendix 0, we take the signal producing region to be O.lSi?©. The Sun's an- 
gular diameter is 32', thus the solid angle VLg subtended by the signal producing region is 



tf and were determined to within 1 s; '60 s' refers to the integration time step for computing the average 
flux. 
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Qs = vr(0.13 ■ 32/2)^ arcmin^ = 13.6 arcmin^. Thus we obtain 

F = 122 A ( _^ _ I . (21) 

\s ^ cm "^keV arcmm / 

Next, we define the unit significance for a single orbit o hke this 

go ^ K^mi^ ^ ^1/2^-1/2 tl/2^$)o ^ (22) 

Q is called quality factor and does not depend on a particular orbit but only on the in- 
strument used for x-ray observation^*^, whereas So is determined by the orbit itself and the 
observation mode. Using Eq. (pTf), we can express Q in terms of /, which yields 

^ , / v4 \ ^/Vs""^ cm^^keV"^ arcmin"2\ , 
Q = 0.09cm^s-/^(^) [ j . (23) 

We now sort all orbits of a particular satellite in decreasing order of Sq and add the first 
/ orbits from the top of the list to the analysis. Thus we can compute the total significance 
S for I orbits: 

S{I)^[Y^{S"f] =[Y.{^oQf] =Q[T.^I] =QS with K = Y.^. (24) 

\i=l / \j=l / \i=l / i=l 

V ' 

=E 

The maximal possible value for I is given by the number of orbits in the covered time period, 
which in our case is one year. The definition of S is inspired by the form of a GauBian y^- 
function and the fact that for 1 degree of freedom, a difference of x corresponds to ^/x a 
significance. In this way, the contribution from a particular instrument, encoded in the 
quality factor Q, and of a particular orbit, encoded in S, can be cleanly separated. The 
only remaining effect of the particular satellite on S is the necessary tcut? where fixed and 
turning mode should bracket most realistic setups. 

For each value of I we obtain a value for S and tu and we can plot S as a function of 
This is shown in Fig. |] for those 3 satellites which have the largest maximal obtainable E 
in turning (red lines) or fixed mode (blue lines). In case all orbits have a very similar value 
of Soi S is approximately proportional to i/t^. The fixed mode satellite 25399 exhibits this 
type of behavior for all times shown in Fig. ^. However, for some satellites (like 13777) all 



-"^^ The background rate F actually has some dependence on the position relative to the Earth and is not 
constant in time. We will comment on this point later in more detail. 
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FIG. 5: Shown is the significance S as defined in Eq. ^ as a function of the observation time 
tu- In blue those three satelhtes are shown for which S reaches the largest possible value in fixed 
mode. Whereas, in red the three satellites are shown for which S reaches the largest possible value 
in turning mode. The big dots denote the times at which each satellite has reached 80% of its 
maximal S in the corresponding mode. The gray shaded area shows S for the case where the 
average flux per orbit is constantly (^)gcx = 4 • 10~^cm^^s^^, which corresponds to the result 
obtained in [Qal. 



available orbits are used up at relatively small values of tu- Obviously, satellites which do not 
have to turn, i.e. the fixed mode ones, have a clear advantage. They can accumulate more 
useful time since they do not lose time by turning once they enter the dark orbit. Moreover, 
as indicated in Fig. ^, the GECOSAX fluxes are highest either at the very end or beginning 
of each dark orbit. Thus, they can reach values of S = 0.0035 cm~^ s~^/^, about 3 times 
larger than turning mode satellites. This factor of 3 is very relevant as 5* only increases as 
the square root of time, area or the inverse background rates, therefore a 3 times larger E 
allows a 9 times smaller area or 9 times larger background while having the same statistical 
significance. 

For turning mode satellites (red lines), there are a few long duration, high flux orbits 
which contribute the bulk to S and afterward the curve increases much more slowly than 
y/t^. Thus for an optimal use of resources, it is advisable to avoid those orbits and to restrict 
the axion search to only the best available orbits. 
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Therefore, we introduce the quantity Ego = O.SSmax, which is just 80% of the maximal 
obtainable S. The reduction in significance is small compared to the savings in observation 
time when the experiment is restricted to those orbits which allow to reach SgQ. Sgo is 
marked by a dot on each curve in Fig. ^ The corresponding times tgo are much shorter than 
the maximal available in case of satellite 29052 the reduction is nearly a factor 10 with a 
minimal sacrifice in E. We see that tgo can be significantly below 10^ s for the best available 
satellites. For example satellite 13777 in fixed mode would reach a Sgo = 0.0028 cm~^ s~^/^ 
in about 1 week^^. The values of Ego, the total time needed and the number of necessary 



orbits in both fixed and turning mode are listed for all considered satellites in Table pi . 

All optimal orbits we identified fulfill the naive expectation stated at the beginning of 
this section: they maximize B^L'^ by having their dark orbits in the regions of the strongest 
geomagnetic field and they all have orbits close to the maximum altitude of 1000 km con- 
sidered here. For the sake of completeness, we have added the performance of the current 
ISS orbit in Table |TTT[ Due to the very low altitude of the ISS of around 380 km, this orbit 
does not perform very well. On the other hand, if the restriction on the maximal altitude is 
relaxed, completely new types of orbits become available. 

One interesting such class are the so called Molniya orbits, which are highly eccentric 
with perigees ~ 1000 km and apogees ~ 40 000 km. These orbits have a period of 12 h and 
have a repeat ground track, i.e. they reach the same point above the Earth every 12 h 
and can thus cross the geomagnetic pole every second orbit. Due to their primary design 
goal of allowing communication with high latitudes in Russia, they have their perigee on 
the Southern hemisphere, some of them very close to geomagnetic South pole. This implies 
that during antarctic night, which corresponds to the summer months on the northern 
hemisphere, a satellite on such an orbit would be in darkness in a very high i?-field region 
at an altitude of ~ 1000 km every 24 h for a duration of about Ih. Therefore, a few of the 



best available such orbits out of 24 tested ones are listed as well in Table |T|. Note that 
these orbits could reach comparable sensitivities to the best available ones discussed so far. 
Especially, in turning mode they could yield a significance 20% better than any other orbit. 
The total observation time needed would be very short, on the order of few 10^ s. However, 
the variability of the Earth magnetic field so far out is greater and would require a more 
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Including the 2 times 60s per orbit. 
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careful consideration than the present note allows for. Also, the obtained bound due to the 
much longer average length of the axion path, would deteriorate at a smaller value of rria 
compared to the calculation presented in figure |. 

VI. SENSITIVITY TO ga^^ 

In order to compute the sensitivity to Qa-y we have to specify a value or range of values 
for the quality factor Q. Clearly, Q is a very instrument-specific quantity and each existing 
x-ray detector in space will have its unique value of Q. However, as shown in Eq. (p2D, 
Q = \J A/F is a combination of two factors: the effective x-ray collecting area and the 
background rate F. In principle, these two factors can be scaled independently. Therefore, 
we will consider the range of effective areas and range of background rates / found in real or 
planned x-ray satellite missions separately. We list the effective area, the background rates 
(/, F) and the resulting quality factor Q in Table |. 



Mission 


ID 


Instrument 


effective area 


background rate / background rate -F 


Q 


reference 








2 

cm 


10 *cm 2 s 1 keV ^ arcmin ^ 


10-'^cm-2 s-i 


cm^ s^/2 




XMM^o 


25989 


EPIC MOS 


900 


29 


36 


5000 


il§ 


XTE 


23757 


PCA 


7000^1 




3600 


2260^2 


nil 


SUZAKU 


28773 


XIS FI 


250 


6.3 


7.6 


5735 




XEUS^o 






50000 


120 


147 


18443 


i3 



TABLE I: Values for effective area A, the background rate / and the integrated 
background rate F for various existing and planned x-ray observatories. For 
those missions already in space we list the US SPACECOM ID number. F is 
integrated over the energy range 1 — 10 keV and the source size of 13.6 arcmin2 . 
Given are also the resulting quality factor Q and the reference for the infor- 
mation. 

Table | is not intended to be an exhaustive survey, but to indicate the possibilities of a few 
contemporary missions. The effective areas range over (250 — 50 000)cm^, the background 
rates F span (8 — 3 600) x 10~^ cm~^ s~^, and the resulting Q values are in the range (2 200 — 

^° Not in low Earth orbit. 

'^^ This number corresponds to the value at the beginning of the mission. 

The PCA is a non-imaging detector and hence there is no background rate / given. Since it covers the 
whole Sun, the axion flux for rs = I has be to taken, which is 1.6 times larger than the one for = 0.13. 
This correction factor has been applied to the Q value quoted here. 
The background cited is the measured value, while observing the dark side of the Earth. 
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18 000) cm^ s^/^. Taking the extreme combinations of the effective areas and background 
rates from this table, the corresponding range of Q is (300 — 81 100) cm^ s^^^. For the 
sensitivity estimate presented in Fig. an effective area of 1 000 cm^ and a background rate 
of 7.6 X 10~^ cm^^ s~^ is assumed, yielding Q = 11 471 cm^ s^^^. The value 7.6 x 10~^ cm"^ s~^ 
corresponds to the background rate measured by SUZAKU while observing the dark side of 
the Earth in the energy range 0.5 — lOkeV. Therefore, this constitutes a guaranteed 



upper bound on any interfering x-ray luminosity from the dark side of the Earth, i. e down 
to this level the dark side of the Earth is certainly dark in x-rays. Note that the instruments 



on board SUZAKU are among the most sensitive ones for extended sources p8 |. 

For the computation of the actual sensitivity to (^a^ we will replace the significance defined 
in Eq. (|2^) by the correct form of x^-function. Since the count rates are very low it is 



necessary to use the Poissonian form of the x^-function, see e.g. [^. There is considerable 
variation in the GECOSAX flux along a single dark orbit as is obvious from Fig. Therefore, 
the signal to noise ratio will also vary greatly and hence this time dependence can be 
exploited. Thus, the data to be fitted consist of the time series of all time bins of 60s 
which are in a dark orbit and are not within t^ut of either tf or t ^ and belong to one of the 
those orbits which comprise SgQ. 

ho,i{Ej) = FAAtAE, 



nZAE,) = KAEj). (25) 

with ^-lo being Qa-^ in units of 10~^° GeV~^ (. . .)mae is the average over the energy interval 
AE and the time interval At. Here, AE = IkeV and At = 60s. Next, we define the 
X^-function as follows 

x\g,o,m^) = 2Y^ ^ T.<L{E„m^-KL{E,) + r^,{E,)\n ^> (26) 

o=l i=l j=l '^thel-'^j) ""^aJ 

The bound gt on gio or ga'y is found by requiring that x^idb, TTT^a) = 4, thus the bound is at 2 cr 
or 95% confidence level. This is repeated for many values of m^. The resulting sensitivities 
for the top performing satellites for both the fixed and the turning mode are shown in Fig. ^ 
as a function of m^. It turns out that accidentally, the same satellite performs best in 
both modes. The asymptotic sensitivities for are 4.7 ■ 10-" GeV"^ for fixed mode 
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FIG. 6: Sensitivity to 5^7 as a function of the axion mass ma at 2 a (95%) confidence level. The 



blue shaded region is excluded by the CAST experiment ||13|] . 



and 6.6 ■ 10 GeV ^ for turning mode. For comparison the CAST asymptotic bound is 
8.8 ■ 10^^^ GeV^^ |13]. Thus a fixed mode observation would improve the CAST limit by 



about a factor of 2, which is considerable given that the signal scales as g^^, i.e. the actual 
performance in terms of flux sensitivity is more than 12 times better than CAST. Given 
the fact that the flux sensitivity scales only as square root of time, area and the inverse 
background rate, the actual increase in time or area would have to be 150-fold to reach this 
sensitivity. 

One important issue for the validity of the result is how it would change if we allow the 
background to have systematic errors and a time variation. From a purely statistical point 
of view, we observe that the number of background events per time and energy bin typically 
is very small with a mean value of 0{1) for the F and A chosen here. Thus the purely 
statistical variation of the background in each single bin is around 100%, i.e. any fully, 
between all time bins, uncorrelated systematic variation of the background would have to 
be of that order of magnitude to produce a visible effect. This seems to be very unrealistic. 
On the other hand, a common mode change of the background, i.e. an effective systematic 
variation of F would be very difficult to accommodate since the data contains bins with 
nearly no signal. Those bins severely constrain F. Therefore, one would have to introduce 
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a systematic error which closely mimics the time dependence of the axion signal. 

Note, however, that almost all time dependence of the background is due to either varia- 
tions in the geomagnetic field, either by position or due e.g. solar wind or due to changes of 
position with respect to the Sun. Thus, typically, backgrounds will be high if the magnetic 
field is low. The signal, however, will be large when the magnetic field is high. This is a very 
strong anti-correlation. Also, the background mostly depends on the total field strength and 
not on the component perpendicular to the axion path. The position of the satellite with re- 
spect to the Sun is not very different for consecutive orbits and hence background events do 
not mimic the predicted time dependence of the signal. The strongest background rejection, 
however, is due to the known direction of the signal. Therefore, the same strategy as used 
by CAST can be applied here as well. The image of the axion-producing core of the Sun in 
the focal plane of the x-ray telescope will cover only a few pixels out of the whole sensor. 
All the other pixels can be used to measure the background in situ. The main systematics 
in that case would be due to pointing errors and the width and shape of the point spread 



function of x-ray optics. In CAST these errors are estimated to be negligible |T3 . 

In comparing the result from a full time and energy binned analysis as defined in 
Eq. ^ and the one of the simplified treatment using the significance S as defined in Eq. |2^, 
we find that using E underestimates the flux sensitivity by about 50%. The resulting error in 
ga-y is about 5%. Thus using HQ to estimate the obtainable sensitivity for seems to 

be a conservative approximation with very reasonable accuracy. Therefore the value of Ego 
in Table fTl| can be used for a prediction of the potential of a given satellite or instrument. 
Using this procedure the sensitivity at A^cx can be obtained by 

* ^ (27) 



10-iOGeV-V V A^ 

Taking the largest Q = 81 100 cm^ s^''^ encountered in the discussion of Table |I| and the 
largest value of S = 3.5 x 10^^ cm^^ s^^^^ found from Fig. |^ we get a hypothetical limiting 
sensitivity of 

gb = 2.9 X 10"^^ GeV"^ (28) 
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VII. DISCUSSION AND OUTLOOK 



The overall accuracy of our GECOSAX flux prediction depends on the accuracy of the 
geometric description, the magnetic fleld model, the air density proflles and their effects on 
x-ray propagation, as well as the numerical implementation of Eq. ([T^). In section we 
concluded that geometry related effects are accounted for within about ±25 km or ±5 s. This 
in turn introduces less than 10% error in the GECOSAX signal prediction. The main sources 
for these errors are, our purely geometric deflnition of the Earth shadow and our treatment 
of TLEs and SGP4. The latter source would be absent in an actual measurement. Also, 
the time of entry into the Earth shadow is actually easily accessible via the telemetry data 
of the satellite, e.g. electricity production in the solar panels should be a precise indicator. 
Concerning the magnetic fleld modeling, we found in section |T| that B-fleld errors should 
be less than 5% resulting in at most a 10% error on the signal. This component may 



be difficult to improve even in a real experiment. In section we found that a static 
average atmospheric model can be safely used without introducing more than 5% error on 
the GECOSAX flux. Numerical integrations and the coordinate transformation should not 
contribute to the total error budget. We verifled our code against available analytical results 



| r?| . We also note that there is an approximately 5% annual modulation of the solar axion 
flux due to the variation in the Sun-Earth distance. This effect was not accounted for in our 
computations and modulates our computed signal at the same level. However, this is not 
a source of uncertainty and is in fact a predicted feature of the signal. Thus, we flnd that 
the results for the GECOSAX flux presented here should have an error not exceeding about 
15%, which in a real experiment may be reduced down to about 10%. Given the above 
considerations, low Earth orbit measurements of GECOSAX provide a novel experimental 
avenue for going beyond the current laboratory bounds on the axion-photon coupling, for 
axion masses below 10~^eV. We hope that the analysis presented here will help motivate 
future experimental efforts in this direction. 
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APPENDIX A: COORDINATE SYSTEMS 



Unfortunately, the different algoritlims and programs used for the calculation use different 
coordinates systems which sometimes have non-trivial transformation properties. Therefore, 
a little digression on commonly used coordinate systems is required. All coordinate systems 
which can serve as quasi- inert ial frames are ultimately defined by astronomical observations. 
The idea is that very distant astronomical objects like quasars allow us to define the orienta- 
tion of a triad in space which does not change with time. This triad then can be attached to 
the barycenter of the Earth. We call this the celestial reference system (CRS). Clearly, the 
CRS is not exactly inertial, but deviations are very small of order 10~* for special relativistic 
corrections and 10"^'' for general relativistic corrections For example, the equations of 
motion of a satellite around the Earth, or of the Earth around the Sun are valid in the CRS. 
Observers typically do not float in space but are attached to the surface of the Earth and 
thus observations will be relative to this surface, which can be used to define the so called 
terrestial reference system (TRS), which due to the rotation of the Earth is clearly not an 
inertial system. 

The task is to find the coordinate transformation from CRS to TRS and its inverse, 
this process is also referred to as coordinate reduction. This problem is, however, very 
much complicated by obsolete notations stemming from times when the distinction between 
astrology and astronomy was not always clear. Moreover, many layers of approximations 
of varying accuracy are present, owing to the difficulty of implementing a full coordinate 
reduction without powerful enough computers. 

Since both the CRS and TRS have the same origin in the barycenter of the Earth^^, the 
full transformation can be described by three time-dependent rotations. In reality, these 
three rotations are often split into four parts, since this makes it easier to derive suitable 
approximations: 

xcRS = P(i)N(t)R(t)W(t)fTRs , (Al) 
where P(t) accounts for the precession and ^(t) for the nutation of the Earth spin axis. 



Note, that the International Celestial Reference Frame (ICRF) has its origin at the barycenter of the 
solar system. All coordinates we use have their origin at the barycenter of the Earth. The transformation 
between our Earth-centered CRF and the ICRF is a simple translation. 
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R(t) describes the rotation of the Earth and W(t) the motion or wobble of the spin axis 
with respect to the surface of the Earth. The precession is caused mainly by the Sun and 
the Moon pulling at the equatorial bulge of the Earth, this is called the luni-solar precession. 
Also, there is some precession of the ecliptic due to the influence of the other planets on 
the Earth orbit around the Sun, called planetary precession. The combined effects move 
the Earth axis by about 50" per year. This effect was already known to the ancient Greeks 
and was supposedly discovered by Hipparchos. Nutation is mainly caused by the fact that 
the Moon's orbit is inclined with respect to the Earth's equator and hence the Moon's pull 
on the equatorial bulges changes throughout each month. The full theory of nutation is 
quite complicated since it receives contributions from many sources. The result is as a main 
period of 18.6 years and an amplitude varying from 9" to 17". 

The Earth rotation is not uniform either and changes e.g. due to the friction caused by 
the tidal bulges. In the simplest case the required rotation angle would be proportional to 

u^t. (A2) 

Instead of defining a time dependent angular velocity, all violations of this simple relation 
are absorbed into the definition of time. The relevant time system is UTl (Universal Time 
1) which is based on observed transit times of distant astronomical objects and basically 
ensures the validity of Eq. |A2| . UTC (Universal Time Coordinated) is the one on which 
our daily life is based on. It is adjusted to keep track of UTl by insertion of leap seconds 
whenever required and never differs by more than ±0.9 s from UTl. 

The wobble of the spin axis exists since the Earth is not a rigid body, but has a liquid 
interior. It is a very difficult effect to predict; fortunately it is a very small effect of only 
0.1". 

In all coordinate transformations used throughout this work, we will neglect effects caused 
by nutation, the difference between UTl and UTC as well as any polar wobble, therefore 
our basic transformation reduces to^^ a^cRS = P(^)R'(^)^trs- This choice of approximations 
is mainly guided by the obtainable accuracy of the satellite orbit prediction system used. 
The errors induced are about ±1 s in timing and less than 30" in angle. 

This choice, in particular, impUes that we do not convert TEME as used by SGP4 into a Mean of Date 
system. 
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Next we need to define the Terrestrial Reference System (TRS) or geodetic coordinates 
(GC). Positions on tlie Eartli are commonly measured by latitude, longitude and height 
above mean sea level. The Earth is not spherically symmetric, but to a very good ap- 
proximation^^ an oblate ellipsoid, with a flattening of about 1/300. This is caused by the 
centrifugal force due to the Earth rotation and the fact the Earth is not a rigid but elastic 
body. There are two consequences from this definition of GC: the vector normal to the Earth 
surface no longer points back to the center of the Earth and the latitude is now defined as the 
angle between the normal to the surface and the equator. The satellite propagation routines 
use the so called World Geodetic System 72 (WGS72), whereas the geomagnetic model is 
based on WGS84. The difference between these two system is less than 6" in latitude, less 



than 1" in longitude and less than 6 m in height |^^. Therefore we use them interchangeably 
in coordinate transformations, whereas the satellite propagation routines use WGS72. 

With the exception of GC, all other coordinates system are just simple Cartesian or polar 
coordinate systems which are related by standard transformations. All numerical algorithms 
follow the description in chapter 3 of . We will therefore give the correspondence of our 
notation with the one of in Table 







Fundamental 


Principal 




Symbol Origin 


Plane 


direction Use 


Notation in Q 


ECI 


Earth 


Earth equator 


Vernal equinox Main system for numerical calculations 


UK 


TRS 


Earth 


Earth equator 


Greenwich meridian Intermediate step in coordinate conversion 


(IJK)iTRF 


GC 


Earth 


Earth equator 


Greenwich meridian Input to geomagnetic model 


LatLon 


TCM 


Site 


Local horizon 


North Output of geomagnetic model 




TC 


Site 


Local horizon 


South Intermediate step in coordinate conversion 


SEZ 



TABLE II: Coordinate systems used and their corresponding names used 
in 1^^. Note that some coordinate systems used are Cartesian whereas others 
are polar. 



We will use Cartesian Earth centered inertial coordinates (ECI) for all numerical calcula- 
tions and thus convert all coordinates into ECI first. The motion of the satellite is directly 



In reality, the shape of the Earth, the geoid, is defined as being an equipotential surface of its gravitational 
potential. There is no simple closed analytic form for the geoid. The deviations from the ellipsoid are 
called undulation of the geoid and are indeed very small < 200 m. Note, that mountains which can be up 
to 10 000 m do not play a role in this paper since x-ray propagation ceases at altitudes well above that, 
thus only the shape of the geopotential iso-surfaces, which determine air density, at heights above about 
50 km are relevant. 
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evaluated in ECI. Also the position of the Sun is directly given in ECI. The Earth magnetic 
field is specified by its location in geodetic coordinates (GC) and the result is a vector in 
topocentric coordinates (TCM). The density of the Earth atmosphere is a function of the 
altitude which is defined in GC. Together with Table |I| we obtain the following chain of 
coordinate transformations for the magnetic field B: 

xeci — ^ a^TRs — ^ -Btcm = B{xgc) B^c — > -Beci , (A3) 

where only the non-standard transformations are given on top of each arrow, which is a 
rotation around the 3 or z-axis to get from TCM to TC. 

APPENDIX B: SIGNAL EXTRACTION REGION AND ANGULAR 

RESOLUTION 

The Sun has an average angular diameter of about 32' and the axion producing region is 
mostly confined to the inner 20-30% of the solar radius. The typical angular resolution of x- 
ray telescopes ranges from arc seconds to a few arc minutes, therefore the Sun is not a point 
source of axions. As a result, the "night-side" image of the Sun in x-rays from GECOSAX 
will cover a finite area in the focal plane of the telescope and it is possible to select a spot 
radius r which optimizes the significance s/\/b. Both the signal and the background will be 
a function of r, which we take be a dimensionless fraction of the solar radius. Since axion 
production is very much concentrated towards the center of the Sun, s steeply rises for small 
values of r and then saturates at r ~ 0.3, whereas the background, assuming it is spatially 
uniform, will rise as r^. Thus there should be a maximum in the significance and the radius 
for which this happens is called r^. Assuming perfect spatial resolution of the telescope and 
no pointing errors, we find = 0.13. In the following we will rescale all values of s/\/b by 
the value obtained at in this case. The rescaled s/\/h as a function of r is shown as the 
black line in Fig. |^. 

So far we have assumed that the telescope has perfect resolution and that there are no 
pointing errors. Both errors have the effect that they will blur the image of the Sun and thus 
the signal density per unit area will decrease, whereas the background density is unaffected. 
Therefore the significance will decrease relative to the ideal case and at the same time r^, 
the optimal signal extraction radius will increase. We assume the point spread function, 
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0.2 0.4 0.6 0.8 1 

r[Ro] 

FIG. 7: Relative significance s/Vb as a function of the solar radius. The different curves are for 
different values of the width of the point spread function as labeled in the plot. 

which contains both the effects from pointing errors and the finite optical resolution of the 
telescope to be a Gaussian with a width or standard deviation of dp. The result of finite 
values for cTp are shown as blue and red lines in Fig. |^. The values of cXp next to each line 
are in arc minutes. The resulting value at maximum is the correction factor which needs 
to be applied to the product in Eq. (^) in order to estimate the limiting sensitivity to 
Qa-y- The SUZAKU telescope has a resolution of better than 2.5' and pointing accuracy of 
better than 0.25' |3^, thus the resulting correction factor is 0.96, which was neglected in 



computing figure 
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APPENDIX C: SATELLITE TLEs AND ORBIT PARAMETERS 



TABLE III: This table lists all satellites used throughout this study. Given 
are the US SPACECOM identification numbers, the name of the satellite and 
its general orbit parameters. The fixed and turning mode columns show 80% 
significance Ego as defined in Eq. the required quality factor to achieve a 
sensitivity as good as the CAST experiment using Sgo> the number of dark 
orbits ngo needed to reach EgOi and the total time needed to achieve Ego (this 
includes the contribution of (cut for each orbit). The last column for each 
mode is the corresponding rank within the satellites in this table. '-' indicates 
that there is no useful dark orbit left after cutting away 600 s. TLE set is the 
number of the TLE set used for the calculations in this paper. All satellite 
data and TLEs are from [M. 



TT~1 
lU 


Name 


Perigee Apogee Inclination 




fixed mode 






turning mode 




TLE set 












Ego Q 

1 


CAST 






rank 


Ego 

CN 

7 


QcAST 


"80 




rank 








km 


km 




1 

a 
u 

1 

o 


CN 
CN 

s 

CO 

o 




CD 

o 




CN 

n 

CJ 
CO 

o 


CN 
Cfi 

B 

CO 

O 




CO 
O 






13777 


IRAS 


888 


924 


80.9 


2.8 


1.2 


588 


0.5 
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404 


20*^22 


COBE 


874 


909 


81.1 


2.5 


1.3 


445 


0.4 


4 












548 




HST 


561 


574 


28.6 


0.7 


5.0 


2301 


4.6 


45 


0.3 


13.2 


1931 


3.9 


35 


70 


^uuoo 


ROS AT 


398 


419 


53.1 


0.6 


5.9 


1754 


3.1 


48 


0.2 


17.1 


1170 


2.2 


38 


293 


21578 


SARA 


725 


754 


81.9 


2.5 


1.3 


2061 


2.9 


5 


0.7 


4.8 


748 


1.2 


3 


148 


21701 


UARS 


368 


475 


57.1 


0.7 


4.8 


1441 


2.4 


44 


0.2 


14.4 


1067 


2.0 


37 


290 


22012 


SAMPEX 


427 


493 


81.6 


0.9 


3.6 


1472 


2.2 


40 


0.3 


11.9 


820 


1.4 


33 


525 


23547 


ORBVIEW 1 (MICROLAB) 


705 


728 


70.1 


1.7 


1.9 


1463 


2.5 


20 


0.6 


5.3 


783 


1.4 


8 


347 


23757 


XTE 


479 


492 


23.1 


0.4 


7.5 


2420 


5.0 


50 


0.2 


20.2 


2090 


4.3 


40 
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25280 


TRACE 


567 
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82.2 


1.9 


1.8 


1418 


1.5 


17 


0.5 


7.2 
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1.2 


17 
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25399 


SAFIR 2 


814 
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81.6 


2.8 


1.2 


1950 


2.7 


1 


0.8 


4.2 
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1.2 


1 
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25560 


SWAS 
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70.0 


1.7 


1.9 


1008 


1.3 


21 


0.5 


7.0 
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1.3 


14 
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25635 


ORSTED 
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83.5 


1.7 


2.0 


2248 


4.6 


24 


0.5 


6.3 


1709 


3.5 


10 
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25636 


SUNSAT 


653 
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83.4 


1.7 


1.9 


2205 


4.4 


19 


0.5 


6.1 
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3.5 


9 
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25646 


WIRE 
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434 


82.7 
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5.2 


47 
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41 
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41 


0.3 


11.1 


1261 


2.4 


30 
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39 
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31 
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25791 
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37 


0.4 


8.2 
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24 
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647 


81.8 
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2201 
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32 


0.4 


7.8 


1714 


3.3 


23 
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25994 


TERRA 


704 


730 


81.8 


1.5 


2.2 


2307 


4.7 


29 


0.4 


7.5 


1854 


3.8 


20 


308 


26033 


ACRIMSAT 


680 


739 


81.9 


1.5 


2.2 


2360 


4.8 


28 


0.5 


7.1 


2003 


4.1 


16 


309 


26546 


MEGS AT- 1 


607 


639 


64.7 


1.5 


2.3 


1384 


2.3 


30 


0.5 


6.9 


884 


1.6 


13 


316 


26561 


HETE-2 


560 


594 


2.0 


0.5 


6.6 


2504 


5.3 


49 


0.2 


17.2 


2219 


4.7 


39 


231 


26702 


ODIN 


582 


611 


82.1 


1.7 


2.0 


614 


0.7 


23 


0.3 


10.1 


285 


0.4 


28 


310 










continued on next 
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TABLE III: continued 
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s 
o 
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26998 


TIMED 


619 


643 


74.1 


1.3 


2.5 


1552 2.6 


33 


0.5 


7.1 


764 


1.4 


15 


943 


27370 


RHESSI 


551 


583 


38.2 


0.8 


4.0 


2049 3.8 


42 


0.3 


11.0 


1638 


3.2 


29 


723 


27598 


FEDSAT 


795 


827 


81.6 


2.1 


1.6 


2179 4.1 


12 


0.7 


4.9 


1687 


3.2 


4 


418 


27599 


WEOS 


792 


827 


81.5 


2.0 


1.7 


2204 4.3 


16 


0.6 


5.3 


1731 


3.4 


7 
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27600 


MICRO LABSAT 


791 
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81.5 


2.0 


1.7 


2206 4.3 


15 


0.6 


5.3 


1730 


3.4 


6 
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27640 
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823 


866 


81.2 
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1.4 
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6 


- 


- 


- 


- 


- 
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27643 
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86.0 


1.2 


2.7 


1359 2.2 


34 


0.4 


9.2 


903 


1.7 


27 
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27651 


SORCE 
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642 


40.1 


1.1 


3.1 


1942 3.5 


35 


0.4 


8.6 


1546 


2.9 


25 


531 


27783 


GALEX 
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29.1 


1.0 


3.4 


2217 4.3 


38 


0.4 


8.8 


1879 


3.7 


26 


331 


27843 


MOST 
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81.2 


2.3 


1.4 
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7 


- 


- 


- 


- 


- 
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27845 


QUAKESAT 


826 
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81.2 


2.3 


1.4 
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9 


- 


- 


- 


- 


- 
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27846 


AAU CUBESAT 
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81.3 


2.3 


1.4 


525 0.4 


8 


- 


- 


- 


- 


- 


163 


27858 


SCISAT 1 


643 


670 


73.9 


2.2 


1.5 


961 1.1 


10 


0.5 


6.7 


913 


1.5 


12 


379 


27945 


KAISTSAT 


677 


715 


81.9 


1.6 


2.1 


2364 4.6 


25 


0.5 


6.5 
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3.9 


11 


347 


28230 


GP-B 
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87.8 


1.8 
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18 


0.4 


7.8 


973 


1.6 


22 
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28368 


DEMETER 
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690 


81.9 


1.4 


2.3 


2301 4.6 


31 


0.4 


7.7 


1822 
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21 
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SWIFT 
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46 
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36 


935 


28773 
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31.5 


0.7 


4.7 
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43 


0.3 


12.3 
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3.8 


34 
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28939 
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13 
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11 
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757 


1.3 


2 
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3.8 


19 
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730 


81.8 
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2417 4.9 


26 


0.5 
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3.8 


18 
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687 
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14 






- 
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29506 
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82.3 
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22 
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The ISS orbit is subject to frequent changes due to maneuvering, therefore this result can only serve as 

a rough indicator of the ISS' orbit quality. 

These orbits exceed the validity range of the geomagnetic model and maybe subject to increased uncer- 
tainties. To minimize this effect, only those parts of the orbit with altitude below IRq are considered. 
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